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ABSTRACT 

We present time series photometry and constraints on additional planets in five of 
the exoplanetary systems studied by the EPOCh (Extrasolar Planet Observation and 



(D ■ Characterization) component of the NASA EPOXI mission: HAT-P-4, TrES-3, TrES- 

^ ■ 2, WASP-3, and HAT-P-7. We conduct a search of the high-precision time series for 



photometric transits of additional planets. We find no candidate transits with sig- 
nificance higher than our detection limit. From Monte Carlo tests of the time series 
p I I using putative periods from 0.5 days to 7 days, we demonstrate the sensitivity to detect 

I Neptune-sized companions around TrES-2, sub-Saturn-sized companions in the HAT- 

' P-4, TrES-3, and WASP-3 systems, and Saturn-sized companions around HAT-P-7. We 

investigate in particular our sensitivity to additional transits in the dynamically favor- 
2 ' able 3:2 and 2:1 exterior resonances with the known exoplanets: if we assume coplanar 

^ . orbits with the known planets, then companions in these resonances with HAT-P-4b, 

WASP-3b, and HAT-P-7b would be expected to transit, and we can set lower limits on 
the radii of companions in these systems. In the nearly grazing exoplanetary systems 
^ I TrES-3 and TrES-2, additional coplanar planets in these resonances are not expected to 

' transit. However, we place lower limits on the radii of companions that would transit 

■ if the orbits were misaligned by 2.0° and 1.4° for TrES-3 and TrES-2, respectively, 

p ■ 
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' — techniques: photometric 
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Introduction 



EPOXI (EPOCH + DIXI) is a NASA D iscovery Program Mission of Opportunity using the 
Deep Impact flyby spacecraft (iBlumd l2005l ). From January through August 2008, the EPOCh 
(Extrasolar Planet Ob servation and Characterization) Science Investigation used the HRI camera 
([Hampton et alJl2005l ) with a broad visible bandpass filter to gather precise, rapid cadence photo- 
metric time series of known transiting exoplanet systems. The majority of these targets were each 
observed nearly continuously for several weeks at a time. In Table 1 we give basic information 
about the seven EPOCh targets and the number of transits of each that EPOCh observed. 

One of the EPOCh science goals is a search for additional planets in these systems. Such 
planets would be revealed either through the variations they induce on the transits of the known 
exoplanet, or directly through the transit of the second planet itself. The search for additiona l 
planets in the EPOCh observations of the M dwarf GJ 436 was presented in lBallard et al.l (|2010l ). 
Because GJ 436 is a nearby M dwarf, it is the only EPOCh target for which we are sensitive 
to planets as small as 1.25 -R©. In this work, we conduct a search for photometric transits of 
additional planets; the tr ansit times of the known exoplanets observed by EPOCh are presented in 
Christiansen et al.l (|201lh . 



The search for transiting planets in the EPOCh light curves is scientifically compelling be- 
cause the discovery of two transiting bodies in the same system permits the direct observation of 
their mutual dynamical interactions. This ena bles constraints on the masses of the two bodies 
independent of any radial velocity measurement (jHolman &: Murravll2005l: lAgol et al.ll2005l ). as has 



been done for the multiple transiting planet system Kepler 9 (jHolman et al 



2oiB). There are also 



separate motivations for searches for additional planets around the EPOCh targets. The search 
for additional transits in the EPOCh observations is complementary to existing constraints on ad- 
ditional planets from photometric observations, radial velocity measurements, and transit timing 
analyses of the known exoplanet. 



Here we briefly summarize such work to date for our five targets. ISmith et al.l (|2009l ) in- 
vestigated 24 light curves of known transiting exoplanets, including the -BPOX/ targets HAT-P-4, 
TrES-2, WASP-3, and HAT-P-7, and found that they were sensitive to additional transits of Saturn- 
sized planets with orbital period s less than 10 days with greater than 50% certainty, although that 



probability is less for HAT- P-4 (IKovacs et al 



timing analyse s of TrES-3b ( O'Donovan et al 



2:1 resonances ([Gibson et al 



because of decreased phase coverage. Transit 



20071 ) have ruled out planet s in interior an d exterior 



20091 ). although the transit times obtained bv lSozzetti et al. (2009 ) for 



TrES- 3b may suggest a deviation from a constant period that could be attribute d to an additional 



body . iFreistetter et al.l ([20091 ) found that a broad range of orbits around TrES-2 ([O'Donovan et al 



20061') would be dy namically stable for additional planets, although the constraints presented by 



Rabusetal 



(120091 1 :or T rES-2 have ruled out an 5 planet in the 2:1 resonance specifically. 



and iHolman et al.l ([20071 ) found no d e viatio ns in the transit timing residuals from the predicted 
ephemeris. Additionally, iRaetz et al.l ([20091 ) observed a candidate transit in their photometry of 
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TrES-2 which might be attributed to a n additional body in the system in an external orbit to TrES- 
2b. However, iKipping &: Bakoa ([20101) investigated the TrES-2 Kepler observations and found no 
unexpected photoinetric decrements and no significant transit timing or transit dur ation variation 



Macieiewski et alj (|2010l ) performed an analysis of the transit times of WASP-3b (jPollacco et al 
20081 ). and found evidence for planet with a mass of 1 5 in a orbit close to a 2:1 res onance with 



the kn own planet. In the HAT-P-7 (jPal et al.l 120081 ) radial velocity measurements, IWinn et al 



((20091) found a drift that provides evidence for a third body. This radial velocity trend is consistent 
with any period longer t han a few months. Finally, the light curves obtained by the Kepler Mis- 
sion (IBorucki et al.ll2010l ) will ultimately enable exquisite constraints on the presence of additional 
planets in two of the systems which were also observed by EPOXI: TrES-2 and HAT-P-7. 

The remainder of this paper is organized as follows. In Section 2, we describe the photometry 
pipeline we created to produce the time series. In Section 3, we describe the search we conduct 
for additional transiting planets. We present a Monte Carlo analysis of the EPOCh observations 
of HAT-P-4, TrES-3, TrES-2, WASP-3, and HAT-P-7, and demonstrate the sensitivity to detect 
additional planets in the Neptune-sized and Saturn-sized radius ranges. In Section 4, we present 
our best candidate transit signals, and from the search for additional transits we place upper limits 
on the radii of additional putative planet in these systems in the exterior 3:2 and 2:1 resonances 
with the known exoplanets. 



2. Observations and Data Reduction 



The photometric pipeline we developed for the EPOCh dat a is discussed at length in | Ballard et al 

(201C) (concernin g GJ 436 in particular), an d is summarized in Christiansen et al. ( 20icl ) (concern- 
ing HAT-P-7) and lchristiansen et al.l (|201lh (concerning HAT-P-4, TrES-3, TrES-2, and WASP-3). 
We outline here the basic steps we undertake to produce the final EPOCh time series. We acquired 
observations of the five EPOCh targets presented here nearly continuously for approximately two- 
week intervals. These intervals were interrupted for several hours at approximately 2-day intervals 
for data downloads. We also obtained for TrES-2, WASP-3, and HAT-P-7 approximately one day of 
"pre-look" observations, implemented to optimize pointing for each target, that predate the contin- 
uous observations by a week. The basic characteristics of the targets and observations are given in 
Tables 1 and 2. Observations of this type were not contemplated during development of the original 
Deep Impact mission; the spacecraft was not designed to maintain very precise pointing over the 
timescale of a transit (Table [2]). Furthermore, the available onboard memory precludes storing the 
requisite number of full-frame images (1024x1024 pixels). Hence, the observing strategy during 
the later observations (TrES-2, WASP-3, and HAT-P-7) used 256x256 sub-array mode for those 
times spanning the transit, and 128x128 otherwise. This strategy assured complete coverage at 
transit, with minimal losses due to pointing jitter exceeding the 128x128 sub-array at other times. 
We elected to exclude the following EPOCh data from the search for additional transits: first, the 
observations of XO-2, for which we gathered only partial transits and had relatively sparse phase 
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coverage due to pointing jitter and data transfer losses. Second, we did not use the observations 
from the second EPOCh campaign for HAT-P-4 (from 29 Jun - Jul 7 2008), which we could not 
calibrate to the same level of precision as the original observations for reasons explained below. Our 
sensitivity to additional transit signals in the HAT-P-4 light curve, which should theoretically have 
improved with additional observations removed in time, was in reality diminished due to the in- 
creased correlated noise in the second campaign HAT-P-4 observations. For this reason, we elected 
to use only the original 22 days of observations in the search for additional transits. 

We used the existing Deep Impact d ata reduction pipeli ne to perform bias and dark subtrac- 
tions, as well as preliminary flat fielding ([Klaasen et al.ll2005l ). We first determined the position of 
the star on the CCD using PSF fitting, by maximizing the goodness-of-fit (with the statistic as 
an estimator) between an image and a model PSF (oversampled by a factor of 100) with variable 
position, additive sky background, and multiplicative brightness scale factor. We then processed 
the images to remove sources of systematic error due to the CCD readout electronics. We first 
scaled down the two central rows by a constant value, then we scaled down the central columns by 
a separate constant value. Finally, in the case of 256x256 images, we scaled the entire image by 
a multiplicative factor to match the 128x128 images (the determination of the optimal correction 
values is performed independently for each target). We performed aperture photometry on the 
corrected images, using an aperture radius of 10 pixels, corresponding to twice the HWHM of the 
PSF. To remove remaining correlated noise due to the interpixel sensitivity variations on the CCD, 
we fit a 2D spline surface to the brightness variations on the array as follows. We randomly drew 
a subset of several thousand out-of-transit and out-of-eclipse points from the light curve (from a 
data set ranging from 11,000 total points in the case of TrES-3 to 20,000 points in the case of 
HAT-P-4), recorded their X and Y positions, and calculated a robust mean of the brightness of 
the 30 nearest spatial neighbors for each selected point. To determine the robust mean, we used 
an iterative sigma-clipping routine that recalculates the mean after excluding outliers further than 
3 sigma from the mean estimate at each iteration (the routine concludes after the iteration when 
no new outliers are identified). Given the set of X and Y positions and the average brightness 
values of the 30 points which lie nearest those positions, we fit a spline surface to the brightness 



Table 1. EPOCh Targets 



Name 


V Magnitude 


Number of Transits Observed'' 


Dates Observed [2008] 


HAT-P-4 


11.22 


10 


Jan 22-Feb 12, Jun 29-Jul 7 


TrES-3 


11.18 


7 


Mar 8-March 10, March 12-Mar 18 


XO-2 


12.40 


3 


Mar 11, Mar 23-Mar 28 


GJ 436 


10.67 


8 


May 5-May 29 


TrES-2 


11.41 


9 


Jun 27-Jun 28, Jul 9-Jul 18, Jul 21-Aug 1 


WASP-3 


10.64 


8 


Jul 18-Jul 19, Aug 1-Aug 9, Aug 11-Aug 17 


HAT-P-7 


10.50 


8 


Aug 9-Aug 10, Aug 18-Aug 31 



^Some transits are partial. 
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variations in X and Y using the IDL routine GRID_TPS. This sphne surface has the same resolution 
as the CCD, and approximates a flat field of the CCD which has been convolved by a smoothing 
kernel with width equal to the average distance required to enclose 30 neighboring points. We then 
corrected each data point individually by bilinearly interpolating on the spline surface to find the 
expected brightness of the star at each X and Y position. We then divide each observation by its 
expected brightness to remove the effects of interpixel sensitivity variations. We used only a small 
fraction of the observations to create the spline surface in order to minimize the potential transit 
signal suppression introduced by flat fielding the data with itself. To produce the final time series, 
we iterated the above steps, fitting for the row and column multiplicative factors, the sub-array 
size scaling factor, and the 2D spline surface that minimized the out-of-transit standard deviation 
of the photometric time series. 

We include two additional steps in the reduction of these data that were not included in the 



Ballard et al.l (j20ld ) reduction of the GJ 436 EPOCh observations. First, during the second cam- 
paign observations of HAT-P-4 and TrES-2, we observed an increase in brightness when the position 
of the star was located in the lower right-hand quadrant of the CCD. At the image level, we ob- 
served a bright striping pattern in this quadrant that caused the measured brightness of the star to 
increase as soon as the PSF entered this quadrant. We found that the dependence of the brightness 
increase in this quadrant was correlated with the Instrument Control Board Temperature value 
recorded for for each image in the FITS header. For the HAT-P-4 second campaign and TrES-2 
observations, we first fit a spline to the dependence of the photometric residuals (after the boot- 
strap flat field was applied, described below) on the Instrument Control Board Temperature, using 
residuals for which the entire PSF of the star fell into the CCD quadrant in question. The most 
egregious brightness increase is 4 mmag, when this effect was most prominent. We then performed 
aperture photometry again on these targets and corrected each image by interpolating the Instru- 
ment Control Board Temperature for that image onto the spline, multiplying this correction value 
by the fraction of the PSF core that fell into the quadrant in question, and dividing this value from 
the photometry. We found that this iterative procedure largely removed this quadrant-dependent 
effect. In the latter half of TrES-2 observations, we no longer observed the brightness increase in 
this CCD quadrant. Therefore, we found that the correction procedure was only necessarily for 
the second campaign IIAT-P-4 and the first portion of TrES-2 observations. However, because of 
the 6-month separation of the second campaign HAT-P-4 observations from the original HAT-P-4 
observations, the behavior of the CCD had sufficiently altered to disallow the combination of the 
data into a single 2D spline surface. The separate 2D spline correction of the second campaign 
observations (spanning only 8 days), coupled with the residual striping artifacts, sufficiently de- 
creased the precision of the second campaign observations that we elected to exclude them from 
the search for additional transiting planets around HAT-P-4. 

Secondly, we include one final correction after we have removed the interpixel brightness vari- 
ations with the 2D spline, which is to perform an additional point-by-point correction to the data 
taken during transit and secondary eclipse of the known exoplanets. The bootstrap flat field ran- 
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domly selects a set of points to create the spline surface, instead of using all the data to create this 
surface; this minimizes the suppression of additional transits. Our sensitivity to additional transits 
is sufficiently diminished during the transits of the known planet that we are concerned more with 
removing correlated noise, and less concerned about avoiding additional transit suppression. The 
two reasons for the diminished sensitivity during transit windows are, first, that we fit a slope with 
time to the points immediately outside of the transit of the known exoplanet (from 3 minutes to 
30 minutes before and after transit) and divide by the slope in order to normalize each transit 



before we fit for the system parameters (this procedure is also detailed in lBallard et al.ll201Cl ). This 
could have the possibility of removing a decrement due to an additional transit. Second, there is 
also the possibility of an occultation of one planet by another. We therefore elected to perform 
a point-by-point correction after the 2D spline for points occurring during the transit and eclipse 
of the known exoplanet, wherein we find a robust mean of the 30 nearest neighbors to each point 
(using the same iterative routine described above) and divide this value individually for each point 
in transit or eclipse. This has the benefit of removing additional correlated noise during transit 
and eclipse, while still minimizing signal suppression of additional putative planets outside of these 
time windows 

After we take these steps to address the systematics associated with the observations, we 
achieve a precision for the unbinned observations which is approximately twice the photon noise 
limited precision for all five targets. We estimate the photon limited precision at the image level, by 
converting the stellar flux from watts per meter squared per steradian per micron (W/m^-sr-/im) 
to electron counts as follows. We first divide by the conversion factor in the FITS header (keyword 
RADCALV = 0.0001175 W/m^-sr-^/m per DN/s), then multiply by the exposure time (INTTIME 



= 50.0005 s), and finally multiply by the gain (28.80 e/DN, per lKlaasen et al.ll2008l for the HRI 
camera). We then estimate the photometric error by calculating l/\/iV, where N is the number of 
electrons. We have excluded read noise, bias, and dark current from the estimation of the photon 
limited precision because these quantities contribute negligibly to the total number of electrons 
measured within the aperture; we briefly summariz e our reasoning here . Using the results of 



the calibration tests on the HRI instrument shown in lKlaasen et al.l (|2008l ). we estimate the read 
noise and dark current (given the CCD temperature of 160 K, as recorded in the image headers) 
to contribute less than a DN. Calculating the median bias value per pixel from the overclocked 
pixels associated with each image, and then multiplying this median value by the number of pixels 
contained within the aperture, we determine that the bias contributes less than 500 DN. When 
compared to the total measured DN flux contained in the aperture, which is of order 10^ DN for 
the dimmest target star, TrES-3, we conclude that read noise, bias, and dark current are negligible. 
We repeat the photon limited precision calculation on 50 images for each target, and take the mean 
of these values to be our estimate for the photon limited precision. Our precision of 1.21 mmag for 
IIAT-P-4 is 94% above the limit, 2.17 mmag for TrES-3 is 106% above the limit, 1.62 mmag for 
TrES-2 is 136% above the limit, 0.97 mmag for WASP-3 is 106% above the limit, and 0.86 mmag 
for IIAT-P-7 is 91% above the limit. The EPOCh precision for GJ 436 of 0.51 mmag was only 
56% above the photon noise limit, which we attribute to the longer baseline of observations with 
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f ewer gaps in phase cove rage, both of which enabled us to create a higher quahty 2D sphne flat 



field (jBallard et al.ll201Cl ). Figure [T] shows five EPOCh time serie s after the 2D spline correction is 



applied; these light curves are identical to the ones presented in [Christiansen et alj (|2010l . |2011| ). 
In the right panel adjacent to each time series, we show how the time series, after the 2D spline is 
applied, bins down as compared to Gaussian noise over timescales of 7 hours (512 points) or less. We 
selected the longest contiguous portion of the lightcurve between transits (and excluding secondary 
eclipse) to calculate the standard deviation as a function of binsize- this unbinned portion typically 
comprises about 2500 points. We compare the expected Gaussian scatter for a bin size of 1 hour 
(assuming that sigma decreases as N~'^^'^, and normalizing to the observed rms of the unbinned 
data) to the measured scatter, and find that the presence of correlated noise inflates the la error 
bar by a factor of 1.86 for HAT-P-4, 1.58 for TrES-3, 1.90 for TrES-2, 2.77 for WASP-3, and 3.14 
for HAT-P-7 for 1 hour timescales. 

We also investigate the transit signal suppression introduced by using a flat field created from 
the out-of-transit and out-of-eclipse data itself. We avoid the suppression of known transits in 
each data set by iteratively excluding those observations (using an ephemeris for the known planet 
derived from the EPOCh observations) from the points used to generate the flat field surface, so 
that we only use the presumably constant out-of-transit and out-of-eclipse observations to sample 
the CCD sensitivity. However, if the transit of an additional planet occurs while the stellar PSF 
is lying on a part of the CCD that is never visited again, the 2D spline algorithm instead models 
the transit as diminished pixel sensitivity in that CCD location. To quantify the suppression of 
additional transits that result from using the 2D spline, we inject transit light curves with periods 
ranging from 0.5 days to 7 days in intervals of 30 minutes in phase (ranging from a phase of zero 
to a phase equal to the period) into the EPOCh light curve just prior to the 2D spline step. After 
performing the 2D spline, we then phase the data at the known injected period and fit for the best 
radius, using as the goodness-of-fit statistic. We show in Figure [2] the radius suppression as a 
function of period for five EPOCh targets. The HAT-P-4 observations occur over a longer duration 
with less gaps in phase coverage, so even at a period of 7 days, we have 95% confidence that the 
radius of an additional transiting planet will not be suppressed to less than 60% its original value. 
For example, an additional 8 i?0 planet orbiting HAT-P-4 will appear no smaller than 0.6x8 i?®, or 
4.8 with 95% confidence. However, for a target with sparser phase coverage, such as WASP-3, 
we have 95% confidence that the radius will not be suppressed to less than 45% its original value. 
The same 8 -Rq planet orbiting WASP-3 will therefore appear no smaller than 0.45x8 i?©, or 3.6 
Rqj, with the same confidence. The average (50% confidence) suppression value of 75% across 
all periods and for all targets refiects the average density of points on the CCD (and thus the 
quality of the 2D spline), which is indicative of the pointing jitter of the instrument. We describe 
our incorporation of signal suppression into our search for additional planets in greater detail in 
Section 3.2. 



-8- 



1.005 
il 1.000 
I 0.995 

0.990 



EPOCh HAT-P-4 Time Series 




Bin Size [Hours] 
0.1 1.0 10.0 



5 10 15 20 160 162 164 166 

Time (Doys ofter First Tronslt Prior to Observotions) 




1.04 




EPOCh TrES- 


-3 Time Series ■ 








i$ 




1 ***** 


! i 




1.02 












1.00 


■ (» 










0.9S 


- ! 


1 


i 




1 i i ■ 



2 4 e 8 10 

Time (Doys ofter First Tronslt Prior to Observotions) 



5 10 15 20 25 30 

Time (Doys ofter First Tronslt Prior to Observotions) 



1.00 „ 



10 100 
Bin Size [Points] 



Bin Size [Hours] 
0.1 1.0 




1.00 „ 



1 10 100 

Bin Size [Points] 



Bin Size [Hours] 
0.1 1.0 10.0 




1.00 „ 



10 100 
Bin Size [Points] 



Bin Size [Hours] 
0.1 1.0 




1.00 „ 



1.00 „ 



5 10 15 

Time (Doys ofter First Tronslt Prior to Observotions) 



10 100 
Bin Size [Points] 



Fig. 1.— Left panels: EPOXI time series for targets HAT-P-4, TrES-3, TrES-2, WASP-3, and 
HAT-P-7. For TrES-3 (second panel from top), we show the hght curve with original modulation 
due to star spots at bottom. Right panels: The standard deviation versus bin size for each target, 
compared to the ideal Gaussian limit (shown with a line, normalized to match the value at N=l). 
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Fig. 2. — The 50% and 95% confidence values for suppression of additional transits as a function 
of orbital period in the EPOCh observations. We have 50% confidence that the transit signal will 
not be suppressed more than the value of the dashed line at that period, and 95% confidence that 
the transit signal will not be suppressed more than the value of the solid line. 
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3. Analysis 
3.1. Search for Additional Transiting Planets 



We search the EPOXI time series for evidence for additional shallow tr ansits. We develope d 
software to search for these additional transits, which is discussed at length in lBallard et al.l (|201Cl ). 
The steps involved in the procedure are summarized in this section. We conduct a Monte Carlo 
analysis to assess how accurately we could recover an injected planetary signal in each of the EPOCh 
light curves. We evaluate our sensitivity to transit signals on a grid in radius and period space 
sampled at regular intervals in B?p and regular frequency spacing in P. We create an optimally 
spaced grid as follows: for the lowest period at each radius, we determine the radii at which to 
evaluate the adjacent periods by solving for the radius at which we achieve equivalent signal-to- 
noise (for t his reason, we expect s ignificance contours to roughly coincide with the grid spacing). 
We use the iMandel &: Agoll (|2002l ) routines for generating limb-darkened light curves given these 
parameters to compute a grid of models corresponding to additional possible planets. If we make 
the simplifying assumptions of negligible limb darkening of the host star, a circular orbit, and an 
orbital inclination angle i of 90°, the set of light curves for additional transiting bodies is a three 
parameter family. These parameters are radius of the planet Rp, orbital period of the planet P, 
an d orbital phase To generate t hese light curves, we also use the stellar radii values determined 
by [Christiansen et al.l (j2010l . 120111 ) from the EPOXI data, with the corresponding stellar masses, 
from the literature, tha t were used to calc ulate those radii. Those radius and m ass values are 



IMRq and 1.26 Mq (|Kovacs et al.1 l2007l) for HAT-P- 4. 0.82 Rq and 0.93 Mq (|Sozzetti et al 



20091) for TrES-3, 0.9 4 Rq and 0.98 Mq (ISozzetti et al.ll2007l l for TrES-2 1.35 Rq and 1.24 Mq 
(|Pollacco et"al]l2008l l for WASP-3, and 1.82 Rq and 1.47 Mq (|Pal et aLlljoOsI ) for HAT-P-7. At 
each test radius and period, we inject planetary signals with 75 randomly assigned phases into 
the residuals of EPOCh light curves with the best transit model divided out, and then attempt 
to recover blindly the injected signal by minimizing the statistic. The period range of injected 
signals is selected for each target individually, to ensure the injected transit signal comprises at 
least two transits in most cases. For a target with high phase coverage, like HAT-P-4, we inject 
signals with periods up to 7 days, but for targets with observations of a shorter duration and sparser 
phase coverage, like TrES-3 or WASP-3, we inject signals up to 3.5 and 2.5 days, respectively. For 
planets with periods higher than these values, we may detect the planet, but with a single transit, 
we expect only a very weak constraint on the period. 

We first conduct a coarse grid search in radius, period, and phase. We select the spacing of 
this grid to minimize processing time while ensuring that the transit was not missed; we polish the 
parameters with a finer search after the initial coarse search. We sample the space at 300 
points in period space (at even frequency intervals between 0.5 and 8.5 days), 50 points in radius 
space (between 0.5 and 5.5 Earth radii) and a variable number of points in phase space set by the 
resolutio n of the transit duration for eac h period. We use an expression for the transit duration r 
given bv lSeager &: Mallen-Qrnelad ([20031): 
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sm I sm 



j — cos^i. (1) 



For each test model, we compute the x^i using the out-of-transit standard deviation to estimate 
the erro r in each point. After the grid minimum is determined, we use the amoeba minimization 



routine ( Nelder Sz Meadlll965l ) to more finely sample the space in order to find the minimum 
from the specified nearest grid point. We also investigate whether aliases of the best-fit period 
from the g^^id improve the fit. We find that roughly half of the best solutions from the grid are 
aliases of the injected period, most at either half or twice the value of the injected period, but we 
test aliases at every integer ratio from 1/35 to 35 times the given period (although aliases other 
than 1:2, 2:1, 3:1, 1:3, 2:3, or 3:2 occurs less than 3% of the time for all targets). We also repeat 
the finer grid search at the three next lowest x^ minima, in case the best solution (or an alias of 
the best solution) lies closer to that grid point. For all injected signals, we recover a solution which 
is a better fit (in the x^ sense) than the injected signal. For this reason, we are confident that we 
are sampling the space sufficiently finely to locate the best solution. 

We quantify the success of this analysis by how well the search blindly recovers the known 
injected transit signal. We define the error on the recovered parameter, for instance period, to be 
I Pinjected — Pobserved \ / Pinjected- Figure [3] shows this relative error in radius, with 95% confidence, 
for all searches. As we note in the last paragraph of Section 2, we anticipate suppression of 
additional transit signals from the bootstrap fiat field treatment of the EPOXI data. We evaluate 
the suppression we expect at the period values used in the Monte Carlo analysis, using the results 
shown in Figure [2l We incorporate this expected suppression by vertically shifting the effective 
radius values of the grid points at which we evaluate our sensitivity to additional transits. For 
example, for HAT-P-4 at 1.63 days, all grid points have been shifted upward in radius by a value 
of 1/0.7, or 1.42, because we anticipate that the radius will be suppressed to no smaller than 70% 
its original value. For this reason, the recovery statistics corresponding to a 3.0 transit depth 
in the final light curve would be accurate for an original transit signal of a 4.3 R^ planet once we 
fold our expectation of signal suppression. 

We also evaluate the overall detection probability for putative transiting planets. Given the 
cadence and coverage of the £^POX/ observations, we determine the number of in-transit points we 
expect for a given radius, period, and phase (where the phase is evaluated from to 1 periods, in 
increments of 30 minutes). We then evaluate the expected significance of the detection, assuming 
a boxcar-shaped transit at the depth of (i?p/i?^)^, and the standard deviation of the time series. 
At each phase and period individually, we scale down Rp to incorporate the signal suppression at 
that ephemeris. We use the improvement in the over the null hypothesis to define a positive 
detection, after we have removed the best candidate transit signal (described in Section 4.1). If we 
do not first remove this signal, then we are a priori defining a "detectable" signal to be any signal 
more prominent than the best candidate signal, and we would be unable to evaluate this signal's 
authenticity. We set our detection limit at an improvement in x^ over the null hypothesis that 
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Fig. 3. — Constraints on radius from the Monte Carlo analysis. For each point in radius and period, 
we create 75 light curves with random orbital phases, inject them into the EPOCh residuals, and 
attempt to recover them blindly. The diamonds indicate the grid of radii and periods at which 
we evaluate our sensitivity; the contours are produced by interpolating between these points. The 
contours indicate the relative error in radius (absolute value of recovered-injected/injected radius) 
that encloses 95% of the results. 
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signifies a correctly recovered period (which we define as a period error of <1%). This number is 
variable among the EPOCh targets due to the precision of the observations and the contamination 
of correlated noise. The detection probability of additional transiting planets, as a function of their 
radius and orbital period, is shown in Figure HI 

For the HAT-P-4, TrES-3, and TrES-2, the A^^ cutoff is set at 250, 200, and 200 respectively. 
For WASP-3, the criterion for detection is 400, and for HAT-P-7, the cutoff is 500. There are 
five exceptions here for these threshold values across the five targets: in our analysis of WASP-3, 
we find two instances of a significance higher than 400, but an incorrect period value: these signals 
both comprise 4 full transit events, and are recovered at a 4:1 alias of the true period of 2.34 days. 
Due to the same instances of correlated noise that produce two positive deviations during two of 
the four transit events that decrease the depths by 2 mmag, a better solution is found at an alias 
of 4:1 that at the true period. For HAT-P-7, we find three similar cases of a 2:1 alias providing 
a better solution than the injected period, although the significance of the detection is above the 
threshold vale of 500. For these three injected signals, which comprise three transits, two of the 
transits are recovered correctly, and the third overlays a single instance of correlated noise that 
decreases the depth of the transit by 0.5 mmag. 

We investigated the true signal-to-noise ratio (including correlated noise) associated with a 
single detectable transi t with the cutoff A y^ significance for each target. We use a method similar 



to the one described by I Winn et al.l (j2008l ) to determine the contribution of correlated noise to the 
standard deviation over a transit duration timescale. We first solved for the transit depth associated 
with the cutoff Ax^ value, assuming a single boxcar transit with standard deviation equal to the 
out-of-transit and out-of-eclipse standard deviation of the unbinned time series. We next found the 
standard deviation at a bin size corresponding to a transit duration for each target. We assume 
an edge-on transit (which assumption we also used for the Monte Carlo analysis) and the shortest 
period where we expect mostly single transits (this period is slightly larger than the largest period 
used for the Monte Carlo analysis; which period range was selected so that we would expect at 
least two transits in nearly all cases). This approximate orbital period is 7.5 days for HAT-P-4, 
4.0 days for TrES-3, 5.0 days for TrES-2, 3.0 days for WASP-3, and 5.0 days for HAT-P-7. Using 
the cutoff transit depth and the standard deviation associated with the transit duration for each 
target, we find that the signal-to-noise ratio associated with the detection criteria is approximately 
constant across the targets, ranging between 5 and 8. The variation in the Ay^ value can be 
attributed in part to the varying presence of correlated noise in the different data sets (and also 
to the number of points associated with each transit, which depends on the transit duration). We 
confirm empirically that planets of these radii are detectable by examining the detection probability 
as a function of radius and orbital period shown in Figure [H We convert the cutoff transit depth to 
a planetary radius, assuming the stellar radius derived from the EPOCh observations and average 
suppression of the radius to 0.75 its original value (roughly constant for all EPOCh targets, as 
shown in Figure [2]). This radius value physically corresponds to the minimum planetary radius 
detectable by EPOCh from a single transit. This value is 7.1 i2© for HAT-P-4, 6.2 for TrES-3, 
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5.3 i?© for TrES-2, 6.5 i?© for WASP-3, and 7.9 i?© for HAT-P-7. Comparing to the nearest radius 
value in Figure HI we find that indeed, at the shortest orbital period where we should expect to see 
single transits, we can detect a planet with radius associated with the detection criteria at high 
significance. At longer orbit periods, we still expect single transits, but the likelihood that the 
single transit occurs during a gap in the phase coverage increases. 



4. Discussion 

4.1. Best Candidate Transit Signals 

We present our best candidate transits here, for each of the five EPOCh targets. Figure [5] 
shows each of the individual candidate transit events that comprise the best candidate signal, as 
well as the entire phased and binned signal. For HAT-P-4, the best candidate is a 2.7 planet 
in a 3.1 day orbit; the Ax^ significance is 61 (as compared to a detection criterion of 250). For 
TrES-3, the best candidate is a 2.9 i?© planet in a 2.63 day orbit; the Ax^ significance is 87 (as 
compared to a detection criterion of 200). For TrES-2, the best candidate is a 3.6 i2© planet with 
a period of 7.22 days; the significance is Ax^ of 269 (as compared to a detection criterion of 200). 
For WASP-3, the best candidate is a 4.2 with a period of 5.9 days; the Ax^ significance is 232 
(as compared to a detection criterion of 400). For HAT-P-7, the best candidate is a 4.4 R^ planet 
with a 3.9 day orbit; the significance of the detection is a Ax"^ of 201 (as compared to a detection 
criterion of 500). The only candidate signal above the Ax'^ detection threshold is the one in the 
TrES-2 light curve; this candidate signal comprises two transit events (the other predicted events 
occur during gaps in the phase coverage) . One of the candidate transit events occurs in a portion 
in the CCD that is never visited afterward; these data are therefore uncalibrated by the 2D spline 
fiat field and are unreliable. Without the transit signal that occurs in the uncalibrated area of 
the CCD, the Ax^ significance of the remaining transit is 80, which is well below the detection 
threshold of 200. 



4.2. Radius constraints 

From the results of our Monte Carlo analysis and phase coverage analysis, we can rule out tran- 
siting planets in the sub-Saturn radius range for HAT-P-4, TrES-3, and WASP-3, the Saturn-sized 
radius range for HAT-P-7, and Neptune-sized radius range for TrES-2. We consider in particular 
our sensitivity to additional planets in the dynamically favorable 3:2 and 2:1 resonance orbits with 
the known exoplanets. In Figure HI we show the detection probability as a function of period for 
planets ranging in size from 3 to 10 R^, with positions of the exterior 3:2 and 2:1 resonances marked 
by vertical lines. We also indicat e in Figur e [H th e regions not guaranteed to be stable by Hill's 



criterion, per the formula given in iGladmanI (|l993l ). Assuming an eccentricity of zero for both the 



known and putative additional planet, the planetary orbits are assumed to be stable if the following 
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condition holds, where /ii = mi/Mstar, = m2/Mstar, a = fii + fi2, and 5 = \J aijay. 



a 



4/3 ■ 



a 



(2) 



We solve numerically for the bound aries in b of the stab le region, using the stellar mas s es and 
masses for t he known planets giv en by iKovacs et al.l (120071') fo r HAT -P-4b, ISozzetti et al.l (120091) 



for Tr ES-3b, ISozzetti et alJ (120071 ) for TrES-2b, iPollacco et all (|2008l ) for WASP-3, and iPal et al 
(j2008l ) for HAT-P-7, and conservatively using a putative mass for the second body equal to the 
mass of Saturn. This results in an overestimate of the extent of the Hill-unstable region for the 
planets with masses smaller than Saturn; while we find we are sensitive to planets with radii well 
below that of Saturn, the mass of putative additional planets depends on their composition and is 
uncertain. However, the critical b values vary slowly with increased mass of the putative additional 
planet, so that increasing the mass to that of Jupiter changes the periods associated with the closest 
Hill-stable orbits by only 7% at most for these systems. In some cases, the 3:2 orbital resonance is 
not guaranteed to be Hill-stable (though it may be stable); the exact boundary of the stable region 
depends on the mass we assume for the additional planet. 

From the detection probabilities shown in Figure HI in the HAT-P-4 system, we are sensitive 
to planets as small as 8 i?© in the 3:2 and 2:1 resonance with HAT-P-4b (with a period of 3.06 
days) with 95% confidence. In the TrES-3 system, with the known exoplanet in a 1.3 day orbit, 
we would have detected a 5 planet in the 3:2 and 1:2 resonance with 70% and 50% probability, 
respectively, and an 8 i?© in either orbit with nearly 100% probability. Around TrES-2 we are 
sensitive to the smallest planets, and would have detected a 4 i?© planet with 65% probability 
in the 3:2 resonance with TrES-2b (which has a period of 2.47 days). In both the 3:2 and 2:1 
resonances, we had a high probability of detecting a 5.0 i?© planet: >95% in the case of the 3:2 
resonance, and 90% in the 2:1 resonance. Around WASP-3, we had 50% chance of detecting a 
5.0 i?0 planet in the 3:2 resonance with WASP-3b (which has a period of 1.85 days), and would 
have seen a planet as small as 8 i?© in either the 3:2 or 2:1 resonance with >95% probability. 
Around HAT-P-7, we would have detected a Saturn-sized 10 R@ planet in either the 3:2 or 2:1 
resonances with 95% probability, and had a 70% chance of detecting an 8 -R® planet. If we assume 
an inclination equal to that of the known exoplanet, we can rule out additional transiting planets 
of HAT-P-4, WASP-3, and HAT-P-7 in the 3:2 and 2:1 resonances of the sizes stated above, as we 
still expect additional planets to transit at those orbital distances. However, the known exoplanets 
in both the TrES-3 and TrES-2 systems are already in grazing orbits, so additional planets in the 
exterior 3:2 and 2:1 resonances would not be expected to transit if they were strictly coplanar with 
the known exoplanet. However, if the orbits of additional planets were misaligned by 2.0° in the 
case of TrES-3 and 1.4° in the case of TrES-2 (using the planetary inclinations and stellar radii 
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then we would observe a transit in both of the 3:2 and 2:1 
exterior resonances. The orbital inclinations of the ice and gas giants in our solar system vary by 
up to nearly 2° (jcoxlboool ). so it is feasible that an additional planet in these systems could transit. 
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Table 2. Characteristics of the EPOCh Observations 



Telescope aperture 30 cm 

Spacecraft memory 300 Mb 

Bandpass 350-1000 nm 

Integration time 50 seconds 

Pointing jitter ± 20 arc-sec per hour 

Dcfocus 4 arc-sec FWHM 

Pixel scale 0.4 arc-sec per pixel 
Sul)aiia\" size 25() X 25() ]>ix('ls si)ainiiii,;^ traiisil. 12^x128 ol liciwise^ 



''With the exception of the HAT-P-4 observations during 2008 January and 
February and TrES-3 observations, which were conducted entirely in 128x128 
subarray mode. 
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Fig. 4. — Detection probability versus period for planets ranging in size from 3 to 10 i?®. The 
detection criteria is set by the percentage of phases at a given period for which the number of 
points observed in transit produces a improvement of the cutoff significance, compared to the 
nuh hypothesis (A^^ of 250 for HAT-P-4, 200 for TrES-3 and TrES-2, 400 for WASP-3, and 500 
for HAT-P-7). We assume a boxcar-shaped transit at the depth of (Rp / R^,)'^ . The vertical lines 
show the positions of the 3:2 and 2:1 resonances with the known planet, and the cros s-hatching 
shows the location of orbits which are not guaranteed to be stable by Hill's criterion per iGladman 
(|l993l ). 
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Fig. 5. — The best candidate transits for the five EPOCh targets. Each of the individual candidate 
transit events comprising the signal are shown at left; the phased and binned signal is shown at 
right. A time of zero on each X axis corresponds to the time of the first transit of the known 
planet observed by EPOCh. The Ax^ significance of the HAT-P-4, TrES-3, WASP-3, and HAT-P- 
7 candidate signals fall below the detection criteria. While the significance of the TrES-2 candidate 
is above the detection criteria, one of the candidate transits (shown in the leftmost panel) occurs in 
a sparsely sampled part of the CCD, so the observations are uncalibrated and unreliable. Excising 
this candidate event, the significance of the remaining signal falls below the detection threshold. 



